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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 1364 

THE PLANE PROBLEM OF THE FLAPPING WING* 1 
By Walter Birribaum 


In connection with my report on the lifting vortex sheet^ which 
forms the essential basis of the following investigations, I shall show 
how the methods developed there sire also suitable for dealing with the 
air forces for a wing with a circulation veiriable with time. I shall, 
in particular, develop the theory of a propulsive wing flapping up and 
down periodically in the manner of a bird's wing. I shall show how the 
lift and its moment result as a function of the flapping motion, what 
thrust is attainable, and how high is the degree of efficiency of this 
flapping propulsion unit if the air friction is disregarded. Finally, 

I shall treat an interesting case of dynamic instability for a spring- 
suspended wing; this phenomenon was confirmed by experiments at the 
Gottingen aerodynamic test laboratory. Professor Prandtl gave me his 
guidance concerning the present report, and I want to express here also 
my sincere gratitude for the abundant stimulation and energetic assist- 
ance he gave to me at all times. 

1. General statements .- The calculations refer to the two-dimensional 
problem, that is, to the wing of infinite length, or bounded by plane side 
walls, or to a wing with so large an aspect ratio that the boundary effect 
is negligible. For the rest, the same assumptions as in the first report 
are valid concerning smallness of the air forces, slight camber of the 
wing, and so forth. The wing is assumed to extend from x = -l/2 to 
x = + 3/2 so that the point x = 0 becomes the center of pressure of 
a plane wing with fixed angle of attack. Let v, in the direction of the 
positive X-axis, be the air velocity at a large distance. Let the posi- 
tive Y-axis point downward. By y = y(x,t), I denote the density of the 
lifting vortices (now a function of the time) . At every variation with 
time of the density of circulation, free vortices of the density e(x,t) 
will separate from the lifting vortices and will drift away with the air 
flow. On the basis of the theorem about the conservation of circulation, 
or by integration of Euler's differential equation once along the pressure 
side and once along the suction side of the wing (compare the unabbreviated 

*"Das ebene Problem des schlagenden Flugels" Zeitschrift fur augewandte 
Mathematik und Mechanik, Band 4, 1924, pp. 277-292. 

■'"Abstract from the Gottinger dissertation of the same title. Avail- 
able in the University library at Gottingen and the state library at 
Berlin. Referent: Professor Prandtl. 

%his periodical (Z.f .a.M.M.) , vol. 3> 1923* PP- 290-297- 
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report) , one finds easily that the vortices must satisfy at every point 
and at every time the following continuity equation of the vortex density^: 


^I + ^£ + v^ = 0 (l) 

dt dt dx 


Generally one will make suitable assumptions regarding 7 (x,t) so that 
e may be found from the equation. If dy/dt is designated by 7(x,t), 
its integral is 


e(x,t) = - 


1 

v 


cp(x-vt) 




"fc + 



( 2 ) 


Where cp is an arbitrary function to be determined by boundary condi- 
tions. After e has thus been found, there results the induced vertical 
velocity in first approximation as principal value of the integral 


w(x,t) = ^ P d | (3) 

2rt J -00 X - £ 


Similarly to the former case, w now has to be put in relation with the 
moved wing contour. If the motion of the wing is indicated by the value 
of the ordinate at the point x, thus by y = y(x,t) (with this formu- 
lation one could also include a moved and simultaneously deforming wing) , 
the kinematic boundary condition reads 

— + v — = w(x,t) (4) 

dt dx 


^Details on the derivation of this equation may be found in a lecture 
of Prandtl "On the Formation of Vortices in the Ideal Fluid" at the hydro- 
aerodynamic conference at Innsbruck 1922. (Lectures concerning the field 
of hydro and aerodynamics, edited by Th. v. Karman and T. Levi-Civite, 
Berlin 192t) and in the original report quoted. 
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This is the equation from which contour and motion of the wing are 
obtained if 7 , and therewith w, is prescribed. In analogy to the 
above said 



If one takes into consideration that the air velocity relative to the 

supporting vortex line has the components v and w - ^ , there result 

at 

as before the air forces for the unit length of the wing 



M 


pv 


1 

'2 


7 X dx 


( 6 a) 


( 6 b) 



The last term in W represents as in my first report the suction 
force at the leading edge of the wing, now of course with a coefficient a 
variable with time of the first fundamental function. 

2. Periodical wing motion .- From these general expressions, I now 
pass to treating the special problem of the flapping wing. It is to be 
expected that 7 itself in this case will be periodical, except for 
a constant mean value. Thus, I determine first from the expression 


7 


r 0 (*)e lvt = roe 1 — = 7 0 e' 


iojvt 


<n'vt 


( 8 ) 
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the free vortices e; for 7 which alone enters into e is independent 
of the constant mean value 7 so that this expression is sufficiently 
general. As customary in oscillation theory, one tacitly deals with 
the imaginary constituent of 7- v is the circular frequency of the 

oscillations, cd = — a dimensionless quantity which I call "reduced 

v 

frequency" (v has the dimensions of an angular velocity since it equals 
the flight velocity, measured with half the wing chord as unit length) . 
If one designates the distance travelled by the wing after a full period 
as the "wave length" X,, one has 


X = 



2n 

CD 


The further results are obtained by series developments with respect to 
the reduced frequency cd which have the form 


2 c mn a;,m ( log u) ) n 


m ^ n 


These series converge satisfactorily only for small cd up to about 
cd = 0.1 corresponding to a wave length of 30 wing chords or more. Thus 
I presuppose for the present report slow quasi-steady oscillations of 
the wing; this does not imply, of course, that the actual frequency can- 
not assume high values, too, if only the wing chord is sufficiently small. 
If I calculate, finally, with cd' = icD, the formulas obtained, at least 
outwardly, real coefficients which offer many advantages for the 
calculation. 


I now assume that parallel flow free from vortices prevails ahead 
of the wing. Then I obtain with x = x - l/2 for e and w, under the 
presupposition that the "switching-on process" of the wing motion has 
already run its course (the motion thus has become steady), the fol- 
lowing expressions: 


e=0; x < -1 


6 = e a = cD'e 


6 = 6n = CD'e^ 


- m .„CD'(vt-x) 


-1 _ _ _ 

* t / _ \ _ . 




6 7 0 U)aft 


-1 £ X < 1 


(vt-x) 


P-1 - 

e a ^7 0 (|)d| 


> ( 9 ) 


x > 1 
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2«w(: 


:,t) . r 2fiUS 
J-l x - i 


t) aT , r +i jr p“ € P (i,t) _ , , 

d S + — — d| + -£ — d| (10) 

j-i x - i Ji x - 5 


As in my first report, I determine 7 Q (x) from the three first 
fundamental functions^ 


7 0 a7 0a + b7 0b + C7 0c 


Where a, b, c now are complex numbers. One has then also 


w = 


V 


.co'vt 


Wo = aw o_ + bw„ + cw. 


Oa 


Ob 


Oc 


The integrals ( 9 ) and (10) cannot be evaluated by elementary functions. 
The last integral (10) in particular leads to an integral logarithm. 

The Euler constant C which appears as a consequence always occurs in 
combination with log 2 so that I shall introduce as transcendent 
number 


Z = log 2 - C = 0.11595 . . . 

especially for the present report. If one takes this fact into con- 
sideration one obtains the w in the form 

w 0a = Oq + + a^oo' log a>' + + a^co'^log a)' + + a£to'3log a)' 

(11) 

^One obtains 


7 0a 


H - x 
1 + x 


7 0b = \I X " x2 7 0c = x fi - x2 
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and corresponding expressions with |3q and are valid for Wq^ 

and w Qc . The coefficients ciq, 3j_, y± are polynomials in x the 

coefficients of which are linear and rational in Z (regarding their 
values compare the original report). Besides, the appearance of log a>' 
shows that a treatment of the problem with a lifting line instead of 
the supporting surface must fail since to the transition to the lifting 
line there corresponds a transition to cjl>' = 0 (v increases arbitrarily 
for decreasing wing chord) . For this reason, the theory of the flapping 
wing here presented is the simplest possible one. 

Bypassing the integral for y, we are now concerned with determining 
the wing motion in an appropriate manner and bringing it into accord with 
the expression for 7. Since the methods of the lifting vortex sheet are 
linear and since a deformation oscillation will not be taken into con- 
sideration, shape and mean angle of attack of the wing may be disregarded: 
It is then sufficient to calculate the oscillations of a plane wing with 
the mean angle of attack 0 and to superimpose these oscillations as 
small fluctuations with time of the profile chord linearly on the arbi- 
trarily prescribed profile. For the amplitudes which are assumed to be 
small, one may disregard the fluctuation of the abscissa of a wing point 
as small of higher order, and has then as the most general motion still 
possible (compare fig. l) 


y = p(t) + xcp(t) (12) 

p and cp are to be developed as Fourier series with respect to ojv; 

I retain only their first term 

y = (A + Bx)^ ,vt (13) 

Higher-harmonic oscillations would again have to be superimposed linearly. 
Fluctuations in the flight direction (which could be expressed by periodica 
fluctuations of v) also will be disregarded. 

I am designating the quantities A and B as complex stroke ampli- 
tude and amplitude of rotary oscillation, respectively, since equation (13) 
is formed by combination of the special translatory and rotary oscillations 

A = 1, B = 0 (a) and A = 0, B = 1 (p) 

Thus 7 q and all quantities linearly connected with it will be linear 

and homogeneous in A and B. For simplification of the calculation, 
the cases (a) and (p) may therefore be calculated separately and after- 
wards be combined linearly, for instance 

a = Aao, + Bap etc. 
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The complex circulation coefficients are found from equation (4) in the 
following manner: a, b, c are expressed in the form of the series 


2 2 

a = a^ + a^jo' + ages' 4°S 0)1 + a, 01 ' + log o>' + 


3 4 


a^ui ' ^log^co ' + agco'3 + ayod'^iog to' + 


agCD'^log^o)' + a^cD'3log3oD' + . . . 


(14) 


b = bQ + . . 


c = Cq + . . . 


Equation (4) has in our case the form 

w(x,t) = (Aco ’ + Bcd'x + B)ve Ui ^ 

= ( aw 0a + bw 0b + CW 0c) eUJ,Vt 


This is to be valid identically in x and t. If both sides are 
interpreted as series in o)' m (log o)') n > their coefficients must be 
identical, a , b, c therein are unknowns. Each coefficient yields 
a relation in x which at first cannot be expected to be identically 
fulfillable by only three free values. It was intended to fulfill the 
condition at least at three suitably selected locations x-j_, X£, and 

Xg. Surprisingly, these conditions are now shown to be precisely of the 
second and not of a higher degree in x. At least this is the case up 
to the terms of third order, and there is no doubt that this will remain 
so for the higher terms as well. The coefficients of second degree 
in x determine, therefore, with their three subcoefficients each one 
triplet of values, each of the a, b, c. In every new comparative 
coefficient there appears a new value triplet of this kind so that the 
thirty unknowns for which a formulation has been set up may be found 
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uniquely from linear relations "by successive evaluation. The calculation 
results in the following values: 


a o 

= 2Bv 



a l 

= 2v j/ 

i +^B(l 

(M 

1 

a 2 

= 2Bv 



a 3 

= 2v£ 

-AZ - B(Z 

- z 2 U 

a^ 

= 2v[/ 

i + B(l - 

2ZJ 

ac 

= 2Bv 




a 6 = 2 v[az2 - b(| + § Z - I Z2 
a? = 2v 02AZ + - Z + 3 Z 2 jj 

a 8 " + " 3Z 3 

a^ = 2Bv 

h 0 = 0 
h-^ = l-Bv 
h 2 = 0 

b 3 = 2 v(a + i b) 
c 0 = 0 

Cj_ = 0 

e 2 = 0 
c^ = Bv 

b.^ = Cj_ = 0, i ^ 4 


(See footnote 5.) 


5 The dimensions of all quantities are affected by the selection 
of the wing chord as unit length. A treatise of the author in the 
Zeitschrift fur Motorluftschiffahrt on the same subject has been arranged 
so that no objections are possible from the viewpoint of similarity 
mechanics . 
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I want to point out here briefly that approximated values for the 
circulation coefficients are obtained also, if the effect of the free 
vortices is disregarded and the elementary calculation made in such a 
manner as if the momentary apparent angle of attack of each wing element 
were decisive for its circulation, namely the angle formed by the air 
velocity v relative to the moved wing element and the direction of the 
latter. Since in case of rotary oscillations (B £ 0) every wing element 
has another vertical velocity, the apparent angles of attack of the 
elements are all different, that is, the wing assumed to be plane behaves 
as if it had an apparent (dynamic) curvature which is periodical. The 
simple calculation (compare the original report, third part, beginning 
of section II) yields accordingly circulation contributions of the two 
first fundamental functions, that is, only the following terms; 



Thus except for higher terms and with consideration of the order of magni- 
tude of B (see below) a good approximation is obtained. 


All the rest follows readily from the circulation coefficients. I 
had introduced the quantity to’ in order to enable an easier calculation 
and consideration also of complex co', that is, damped or excited oscil- 
lations. It is true that for damped oscillations the integrals (10) lose 
their meaning, since the oscillation had been assumed to have been going 
on for an infinitely long time so that here infinitely large amplitudes 
would have had to precede. However, the formulas obtained may be retained 
as approximations, provided a correction is made for the starting process. 
For the following consideration of the steady oscillations, I continue the 
calculation for reasons of physical clarity with the real frequency o>. 


Lift and moment are purely periodical quantities, that is, their 
temporal mean value is 0. Nevertheless, their amplitude and their phase 
angle are of interest. According to formula (6), there is 


K = pv 2 Jt^Ak a + Bkp)e iajv_t = pv 2 Jt^Ak a + Ck^e^ 07 ^ 

= pv 2 7ta>^Ak a + Cky) 

M = pv 2 it^Am a + Bmp) e iiwt = pv 2 it(Am a + Ciay)e iaJV_t 

= pv 2 JttD (Ai^ + Cm^) e iawt 


( 15 ) 
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Therein 


B = Coo; = ockp; 


1% = <“%,» k 3 


I have introduced here the quantity 
lation. This was done because in case of ordinary flapping motions 
B is of the order of magnitude oA as a simple consideration shows. 

The value of o> is assumed to be small; thus A and C will be of 
the same order of magnitude so that calculation with C instead of B 
will be more convenient in practice. The third form of the air forces 

permits, for cases of equal stroke velocity Aiae^ 07 ^, comparison of 
these forces in a simple manner; the roughest approximation theory would 
yield constant coefficients ^ = 2i, k^ = 2. The coefficients k and 
m are complex numbers the constituents of which are given by the fol- 
lowing series developments . It is noteworthy that the series for m 
are finite. 


xay = amp; = a 4^; 

= mp = 


C as new amplitude of rotary oscil- 


k 0 = V + *V; “a = V + V» a = 


= k^ ' = -(l - 2Z)ao 2 - 2a^log a> - 2rtZuj3 + 2itaa3log a> + . . . 


a*a" = V = 2cn - 


jtar 


+ — ZT^j a,3 + 4Za)^li 


og a) - 


> (16) 


2cD^lOg 2 0) + . . 
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tp ' = ky ' = ± ky ' = 2 - no) + {$L + 2Z - 2Z £ 


"*) 


co 2 - 


2(l - 2Z)co 2 log a> - 2co 2 log 2 co - (— - £ + nZ - 3itZ 2 jco3 + 


n(l - 6z)a^log co + 3>Ki>^log 2 co + 


3" = ky" = i ky" = (3 - 2Z)a> + 2co log co - it(l - 2Z)co 2 - 


> (IT) 


2nxrlog a) + 


(* 2 + x 

It + ? 


g + Z - | Jt 2 z - Z 2 + 2Z 3 )co 3 + 


^ n 2 - 1 + 2Z - 6z 2 ^co 3 log co - (l - 6z)co 3 log 2 co - 


2a>^log3i 


'co + 


ma* = ™ia' = - \ aa 2 %" = 0 ®^" = 0 (l8) 


"V = m y' = ffi V = “ | “ 2 mp” = my" = i my" = co (19) 

It suggests itself to represent kg- and nig by the initial posi- 
tion of "time vectors," visualized as rotating, in the complex number 
plane (as customary in alternating-current techniques) and to combine 
from them - with the parameters A and C (of which A may be assumed 
real without impairing the generality) - linearily, in the known manner, 
the amplitude coefficients k = Al^ + Cky (and correspondingly m) 

according to magnitude and phase. The diagram (figs. 2 and 3) shows the 
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curves of the end points of the vectors k and m as functions of the 
parameter cd. In the representation of the curves for k and m which 
would yield the most accurate values for the graphical evaluation, the 
curves of the various approximations have been plotted side by side for 
comparison of the convergence of the series . Of particular interest is 
the case where the wing - without being affected by significant air 
forces - glides over an undulating streamline course, clinging to it 

as much as possible. To this corresponds the parameter C = -1A, which 

in fact yields small air forces of second order in o>, namely 

i k' = 2 cd 2 + . . . | k" = (i - 2zja> 3 + . . . 


1 

A 


m T 


— u? 
2 



3 . The induced drag .- The induced drag is no longer linear in the 
circulation so that the complex method could not he retained without 
new stipulations. It offers no longer any simplifications, and it is 
advisable to continue from here on the calculation with the imaginary 
constituent of all quantities in real form. If the calculation is 
carried out according to equation (7)> W assumes the form 

W = Wq + W-j_ sin(covt + cp]_) + W 2 sin(2ajvt + 93 ) 

Here Wi is different from zero only when the oscillation is superimposed 
on a const ant angle of attack different from zero. W is in oj small, 
of second order. The purely periodical terms are, therefore, hardly sig- 
nificant; however, the temporal mean value Wq, which is different from 
zero is important. I calculate only this value and write for it, for 
reasons of simplicity, again W. I equate 

A = A' + 1A", B = B' + iB", C = C’ + iC" 


and may assume, without restricting the generality. A" = 0. Then W 
becomes a quadratic form in A' = A, B’, B" (or A, C*, C"): a simple 

deliberation shows that the coefficients of B'2 and B"2 (or C'2 
and C" 2 ) are equal and that the coefficient of B'B" (or C'C") is 
zero. Thus W becomes 


W = pv^jA^ + 2AB-w ap ' + 2AB"v aP " + (B’2 + B" 2 )w p p} 
= pv 2 rt {A^ + 2 AC’w ar ' + 2AC"w ar " + (C 2 + C"2)w 77 J. 




( 20 ) 
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0 


The w ik 8X6 a 6 ain series in cu^co log a)) m and have up to higher terms 
the values 


v aa ~ -a)2 + ” (4 " 2ZhAlog o> + o>^log 2 cD + . . 

w a3' = ~ v a,y' - “ £ (3 - 2Z)cd 2 - | a> 2 log a> + 

i (3 - 2Z)a? + $ a^log cd + . . . 

4 2 

V ■ 33 w ar" - - I <“ + f »“ 2 - (I * 2 - \ z2 + I Z P + 

i (l - 2Z)a>3log a) = i a^log 2 <D + . . . 

W 3P = ~2 w 77 = \ ^ ■ \ ^ + 1 ) a)2 + g (3* 2 + 2 - 4z 2 - 4z)cu3 

— (l + 2Z)a^log 00 - — a)^log 2 <s + . . . 

2 2 


> ( 21 ) 


W may be positive or negative, according to the selection of parameters. 
Depending on the type of motion, one has, therefore, to expect drag or 
thrust. I postpone detailed discussion until after calculation of the 
power requirement and the efficiency. The case of gliding over an 
undulating flow mentioned above, corresponding to C ' = 0, C" = -A, 
gives 


W = 


.2 2 

A pv 


n -a 


-oa^log 


CD 



> 0 


Thus the selection of parameters made does not yet correspond exactly to 
the case W = 0; a small correction would have to be provided for this 
case. 
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4. Work done at the wing .- In the free vortices behind the wing, 
energy is contained which must be produced by mechanical work on the 
airplane. This may be done in two ways. The flapping motion may result, 
as mentioned before, in positive or negative thrust. In the case of neg- 
ative thrust a propeller which overcomes this and all other resistances 
to flight is required for maintenance of equilibrium of motion. In the 
case of positive thrust a propeller is needed only until the thrust due 
to flapping exceeds the resistances to flight, whether the flight be 
uniform or accelerated. As to the work performed at the wing itself, the 
wing motion consumes, of couse, energy if thrust exists; the ratio of 
thrust power -L^ = -Wv and the total mechanical power Lf to be applied 
to the wing may then be denoted as aerodynamic efficiency of the flapping 
wing. In case of drag, two more possibilities exist. First, the flapping 
motion may require additional work . The efficiency defined above then 
becomes negative and arbitrarily large when the wing power Lf decreases 
more and more. Second, the case may occur that Lf becomes negative, 
that is, the wing then is supplied with energy from the air (indirectly 
by the propeller) and may use that energy for surmounting the resistances 
in the oscillation mechanism, or may store it in the oscillation itself, 
that is, increase its amplitudes. Aside from this "incremental power" 
the propeller must, of course, in this case yield additionally the energy 
of the free vortices so that one may define the quotient of -Lf and 
the total propeller power as efficiency referred to the power absorption 
of the wing. This is then exactly the reciprocal value of the efficiency 
defined above which in this case, as the quotient of two negative numbers 
becomes positive but larger than one, thus loses its physical meaning. 

Lf is divided into two parts . One has 


— = p + cpx = i<nv(A + Bx)e icDV ^ 
dt 


Therewith Lf = Kp becomes the "flapping power" and L r = Mcp the power 
opposed to the rotary oscillation; thus the wing power Is Lf = Lf + 1^. 
The power opposed to the drag is denoted by Ly- = Wv. In the air there 
then remains in all L=Lf+L w = Lf+L r + L w . I indicate of all Lf 

again only the temporal mean values for which I obtain quadratic forms 
of the same type as for W 


Lf = PV 3 * |A 2 Z iaa + 2AB'* ia p. + 2AB"Z lap ,. + (B« 2 + B"2) l ±pp J- 
= pv3* |A 2 Z iaa + 2AC'Z ia7 , + 2AC"l iar , + (C'2 + C"2) 


V (22) 
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L r has finite series and is, in general, small compared to Lj.. The 
coefficients are individually 

Z vik = w ik (23) 


toux 


-co - 2 + " 2^co^ + 2Zco^log co - co^log^co + 


Z tap' = J z ta7' = J (3 - 2Z)co 2 + | co 2 log 


CO - 


^ (1 - 2Z)co^ - co^log co + . . . 


J tap” = 


1 W* - h “ f ^ + (t + 1 z ■ I z2 


CO 


IP 


i (l - 2Z)a3log co - i co^log 2 co + . . . 


> (24) 


Zt 33 ^2 H 77 - 0 


^raa “ 0 


Z mp' “ “ l ray' ~ 0 


Z rap” = “ Z ra7" = | 


( 25 ) 


1 2 

fc rpp “ ^2 ^77 = 2 “ 


= — ?r 
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haa. - “ 2 - f “ 3 + n, 


— - Z 2 )cjd^ + 2Zo/ 4 'log <D - cn^log^co + 
v 4 / 


l fap' - 5 ! far’ - 5 <3 - 2 Z)a > 2 + | o> 2 l°g a 


A (l - 2 Z)o 3 - A a^log a) + . 

4 2 


W . i W = 5 “ - J “ 2 + (t + I z ■ I z2 )“ 3 


1 (l - 2Z)cD^l0g 0) -i co3l0g 2 0) + . . . 

2 2 


, 1 , 12 
Z fpp = — l *77 = g “ 


CD 


it 2 « 4 

l__ = — or cd + . . 

aa 2 2 


W - | V' - f “ 3 - 


. 1 , it 2 it* /1V 3 . 

V’ - - V = 2 ® " T + 


pp 


= i. Z = - CD - — CD 2 + *3 + s. _ 7L Z 2 - A z\cd^ + 

^2 77 2 2 \^8 4 2 2 / 


— (l + 2Z)cD^iOg CD - — CD^log 2 CD + . 

2 2 


> (26) 


( 27 ) 
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L, as the energy of the vortex trail, can of course never become nega- 
tive and must therefore be a positive quadratic form - definitive in the 
amplitudes. The fact that L in the form here noted is capable also of 
small negative values is caused by the neglect of higher terms in the 
series development. What is obtained in this case, is therefore only the 
error, accidentally negative, of the almost vanishing vortex power. 

Since, when C is used, all coefficients contain the factor a ? , the 
latter has been cancelled in the graphical plotting which is in agree- 
ment with the presentation of the power for constant flapping velocity 
(fig. h) . The curves show that the essential terms always stem from the 
stroke amplitude A, possibly in combination of the latter with the 
amplitude of rotary oscillation, and that the corresponding coefficients 
increase somewhat more slowly than 0 ? . By rotary oscillation I meant 
above an oscillation about the z-axis. More generally; every oscillation 
where the ratio of A and C is real is a rotary oscillation about the 
fixed axis with the abscissa a where A then is A = -aCto. It is 
shown that for not too large values of aoj, that is, for axes which do 
not lie at too great a distance, W and L axe always positive; that 
is, it is not possible to obtain thrust or power absorption by rotary 
oscillations about fixed axes. Production of thrust always requires a 
stroke amplitude different from zero, power absorption requires addi- 
tionally a rotary oscillation lagging by about 90°. Pure stroke oscil- 
lation without rotation also produces thrust which in roughest approxi- 
mation results as -W = , similar to the so-called Knoller-Betz 

effect (if one calculates with the y-axis as the "polar" in the plane 
problem) . 

One obtains good insight into the variation of drag and power if 
one varies, for fixed absolute value of the amplitude ratio C/A = c, 
only the phase angle cp between the two oscillation components where 
one then has to put 


C' = Ac’ = A | c J cos cp, C" = Ac" = A J c | sin cp 


W is shown to become a minimum, the thrust thus a maximum, when the 

w " 

rotation leads by somewhat more than 90°, namely by <p = arc tan — — . 

*0,7’ 

It is plausible physically, too, that the thrust will assume large values 
precisely then when the phase is shifted by about l80° compared to the 
phase which is present for gliding free from air forces over the wave 
course. For every |c| there exists an cd and vice versa for which 
the thrust maximum is absolute. One then has 
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c ' (cd) 


v ay 1 
W 77 


c"(u>) 


w 


ay 


w. 


77 


ft 


^min _ A^pv^Trw m j_ n 


w min “ w aa “ 


v ar ' 2 + ^7” 2 


w. 


77 



L becomes a minimum and disappears for suitable cd except for terms 
of the fifth order when the rotation is lagging by somewhat more than 90 °. 
Finally, the wing power becomes a minimum - thus, the absorbed power 

a maximum - when the rotation lags by somewhat less than 90 °, namely for 


cp = 


arc tan 


^fqy" 
z fa 7 ' 


This maximum too becomes absolute when between 


and cd the following condition is satisfied 


c 


c'(cd) = - 


c"(cd) = - 


^fay ' 

l ?77 

l fay" 

l f yy 


L fmin - A 2 pv5jtl fmin 


^fmin ~ l faa 


l, fay' 


+ l 


L fyy 


fay" * 

= -1 + i CD + . 

C 


Altogether, L becomes negative only when |c| > 1. 
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The efficiency 



Lw 

3=1 

— — where 


is valid for 

L f 

Lf 

^2 

for L f < 0. 



(28) 


is the quotient of two quadratic forms and capable of a great many values. 
Generally, it becomes negative for small or completely vanishing stroke 
amplitudes, and arbitrarily large with co — ^ 0. The same is valid for 
rotations about a fixed not too remote axis. If, however, A £ 0 and 
also lim A £ 0, tj approaches for o> — >0 to the limit 1. 
o>=0 

By a rotary component (c” = 0) of equal or opposite phase tj is 
always deteriorated compared to c = 0; the same is true for c' = 0, 
c" > 0. In contrast, the efficiency is improved for c ' = 0, 

-1 S c" S 0, as can be seen from the diagram (fig. 5). For c' = 0, 
c" = - 1 , tj becomes identically 1, and for c' =0, c" < -1 there 
results power absorption. The representation for fixed jcj in depend- 
ence on cp as a discontinuous single-wave-harmonic function is very 
graphical for the efficiency as well. Figure 6 shows clearly at what 
phase angles the transition from power absorption to power production 
takes place. The most important ones among the coefficients found from 
the series developments have been compiled in the numerical table 
(table I). 

5. Application to the flutter of elastically supported wings .- The 
derived laws could be practically applied in the investigation of a 
phenomenon our pilots observed in the last war. In the so-called sesqui- 
planes, the lower wing was fastened to one single spar only, thus was 
only slightly elastic against small deflections and rotations. In case 
of increased flight velocity, for instance in steep dives, there occurred 
sometimes vigorous flutter of the lower wing tips which underwent obviously 
unstable oscillations in the increased air flow. Of course, such unstable 
oscillations are possible only if energy is supplied to the oscillating 
system, and this occurs, according to my investigations, only when the 
vector of the amplitude of rotation lags by about 90° and when the ampli- 
tude of rotation itself is sufficiently large ^|c| must be > |A Q . Let 
us visualize again the gliding - almost free from air forces - of the wing 
over an undulating flow course where the airspeed (relative to the wing 
elements) has no vertical component. If the amplitude of rotation is 
smaller than corresponds to this case, the motion is damped by the 
co vinter acting air force. If the amplitude of rotation is, on the con- 
trary, larger than in the case above and the wing therefore scoops more 
deeply into the air, the air force always acts in the direction of the 
motion, and the motion is excited. 
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I consider a wing supported on spars, elastic with respect to trans- 
lation in the y-direction and with respect to twist . In order to be able 
to go on from my formulas used so far, I introduce the directional forces 
per unit length in z-direction; I calculate therefore as if these forces 
were distributed continuously over the length of the wing. This assump- 
tion does not lead to any contradictions if the wing in itself, aside 
from its support, is sufficiently stiff. Since I disregard deflections 
in the flight direction, I can show that a wing supported on spars always 
has only three essential elasticity parameters, corresponding to the 
three constants of the work of deformation quadratic in p and cp. 
(Compare fig. 1 and the original report.) With respect to its elastic 
properties, this wing may therefore always be replaced by a wing which 
is supported only on one spar (the "elastic axis") with the abscissa a, 
elastic with respect to translation by means of the directional force c, 
and with respect to rotation by means of the directional moment 7, as 
schematically indicated in figure 1 . The resultant of the elastic forces 
and its moment at the origin then are 


K = -c(p + acp) 

M = -cap - (ca 2 +7)9 

If s is the abscissa of the center-of-gravity axis (point S, 
fig. l), the momentum and its moment at the origin have the values (with 
the mass m per unit length and the corresponding moment of inertia 
■d = mr 2 for the center-of-gravity axis) 

G = m(p + scp) 

U = msp + (ms 2 + d)cp = msp + m(s 2 + r 2 )cp 

The equations of motion of the wing then read 

G = K U = M 

For K and M the air forces have to be added. If I put as before 
p = Ae^ cp = Be“ ,vt 
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and take into consideration that the air forces have on the wing the 
opposite effect from the one the forces of the wing calculated above 
have on the air, K and M become 


K 


_ _ e a>'vt -jcA + caB + pv^^A!^ + Bkp)|' 

M = -jcaA + (ca 2 + y)B + pvMAmx + BmpH 

For abbreviation, I introduce the following designations 



c 

mv 





Therewith there result finally as the 


equations of oscillation 


A(o>' 2 + o^, 2 + Pok a ) + b(scjd' 2 + aac^ 2 + p Q kp) 

A(sa>' 2 + aaiQ 2 + pQm^) + B^(r 2 + s 2 )u)' 2 + (a 2 + q 2 )o^ 2 + 

Equating the determinant to zero yields the equation for the fre- 
quencies a) = -ico'. If I first disregard the forces (for this purpose 
I put Pq = 0) , there result two main oscillations as a consequence of 

the coupling of the oscillation of the center of gravity with the fre- 
quency cjjq and the oscillation about the center of gravity with the 

frequency ^ oDq. With 5 = (s - a) 2 + r 2 + q 2 there is 

u> 1>2 2 = -u 1>2 ' 2 = “\( & ± 'Js 2 + ^i 2r2 ) 

2r 


= 0 

pQmp) =0 (29) 


A 

B 


s(d t ]j 5 2 - 4q 2 r 2 j - 2ar 2 
5 ± \ 8 2 - ^q 2 r 2 - 2r 2 




( 30 ) 


y k = B^x - a^e 1 ^^; k = 1.2 
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The main oscillations therefore are rotary oscillations about fixed 
axes at the distances a^a^. The quantity s - a forms a measure for 

the coupling. In general, there develop beats from both main oscilla- 
tions. If I write those in the form 


y = 


,a l B l e 


ievt 


- a 2 B 2 e 


-ievt 


+ x 


[Bx‘ 


ievt 


Bg e 


-tevt 


3 } 


,ioovt 


oq + cr>2 = 2od 


o>p - erg = 2e « a) 


The motion may be interpreted as an ordinary oscillation with an ampli- 
tude ratio slowly variable as to magnitude and phase. Therein there 
appears of course, periodically recurrent, the phase angle which corre- 
sponds in the air flow to the power absorption. If the air forces are 
to counteract the change of this phase angle , corresponding to con- 
tinuous supply of energy, it can be shown that the case of slight coupling 
(s - a small) for balanced or almost balanced frequencies of the uncoupled 
system is the best presupposition for this phenomenon. I am anticipating 
from the results of the following calculation with consideration of the 
air forces that without coupling (s = a) no unstable oscillation at all 
would be possible. 

The air-force coefficients occurring in the oscillation determinant 
are themselves functions of the frequency 03 f which cannot be indicated 
by simple analytical expressions so that the roots of the oscillation 
determinant cannot be obtained in a simple manner. However, since their 
existence is secured by the physical meaning of the problem, it is per- 
missible to introduce into the equation instead of the coefficients k 
and m, the first terms of their series developments, all the more so 
since the occurring factor pq generally is a small number; for it is 

sensible to break off the series for cos x in order to find from the 
polynomial obtained for instance approximately the first zero of this 
function whereas the same method is meaningless for e x , since no zeros 
of this function exist in the finite domain. 

The oscillation determinant obtains the following form 
r 2 o) f ^ + 5 <Uq 2 cd i2 + q 2 a^ + P 0 a) ,2 j7^r 2 + s 2 )!^ ” sk (3 " sm a + m (T] + 
p o“o 2 H a2 + <l 2 ) k a - + “g] + Po^karap “ = 0 


(3D 
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The roughest approximation is the result of the calculation which 
bases the determination of the lift and moment coefficients - in the 
manner of the theory of the Knoller-Betz effect - on the momentary 
apparent angle of attack and the apparent (dynamic) curvature of the 
wing. Even with this procedure, there result in certain cases complex 
roots cd' with positive real parts which result in an "increment" of 
the oscillations and thus correspond to dynamic instability. I shall 
not here discuss this approximation more closely. Also, I shall only 
briefly mention the case of the greatest instability since this case 
is trivial. It occurs if the elastic axis lies so far to the rear that 
at the slightest displacement of the wing from equilibrium position the 
air flow simply causes the apparatus to tip over aperiodically toward 
the rear. This always occurs as soon as a 

velocity sufficiently large, namely for a > 

for a > 0. 

If I retain of k and m all terms up 
the period-equation of the wing becomes 


is positive and the air 

< i 2a \) 2 _ 2 . y 


or v^ > 


2P, 


0 


2pjta 


the second order in a>'. 


f(co’) = 


With the 


+ aQ'cn'^log <d' + aQ^'^log^jb' 
CqCJo'^ + CQ'cn'^log co* + Cq " 00 logics' 
dQ'ui' log co' + e^ = 0 
coefficients 


+ bQ0)'3 + bQ'o)'3lpg co' + 

+ IqCu ' + 

(32) 


a 0 


= r2 * p 0 £ . p 0 2 (i 

b Q = P 0 P + P 0 2 (i - z) > 0 


> 0 


c 0 " “0 (p + P 0 V ) “ 2p 0 S + P Q 2 


d o “ ^0 p 0^ > 0 
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e o = “o 2 ^“o 2 “ 2p O a ) 

a 0 « = P 0 [2(r 2 + s 2 ) - 2s(l - 2Z) - p Q (± - 2z) 

a 0 " = -p 0 (2s + p o ) 

V = " p o ( 2s " p o) 

c 0 ’ = EpqO^P + q 2 - a(l - 2Z)J 

c 0 ’’ = <V = - 2 P 0 au \) 2 

S=(s-a) 2 +r 2 +q 2 >0 
(i=l + 2(r 2 + s 2 ) _ s(3 - 2Z) > 0 

\ = 1 + 2(a 2 + q 2 ) - a(3 - 2Z) > 0 

e = ( r 2 + s 2)(l - 2Z) - s(l - 2Z + 2Z 2 ) + | > 0 

v = (q2 + a 2 ) (l - 2Z) - a(l - 2Z + 2Z 2 ) + | > 0 

i 


Z = 0.11593 . - 
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6. Numerical evaluation . - The equation can be solved only approxi- 
mately, of course. For this purpose, I first omitted the logarithmic 
terms and determined the roots do' of the algebraic equation 


g(oi') = agCD' 1 *' + bgCJO'^ + C ^' 2 + dgCO' + e Q = 0 (33) 


I regard this equation as an approximation, equate 03' =03' + x, and 
can now develop the logarithmic terms retaining the linear terms in X- 
Thus I obtain 


X = -03' log 03' | ( 34 ) 

R = a 0 'a3' 3 + a 0 "a3'3iog 03' + b^d? 2 + Cq'cd' + d Q '(l + 03' log a?) 

S = 4agd3'3+ aQ'a3'5(i + 4 log 03 ') + 2a 0 "d3'3iog oo'(l + 2 log go') + 

3bQ03' 2 + b 0 'd3 ,2 (l + 3 log do) + 2 cq(U3 ’ + Cq'o3'(1 + 2 log So') + 

d 0 + cLq'(1 + SdS' log do')(l + log do') 

The procedure may be continued and yields the further approximation 

t R(03')a3' log 03' + g(o3 ' ) 

X' = £(^j 

Finally, there follows the complex amplitude ratio B:A = b from 
one of the equations ( 29) . 
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Finally one now has to learn the conditions for which the equation 
for a)' has complex roots with positive real part, corresponding to 
excited oscillations or critical support of the wing. For the limiting 
case of dynamic indifference, I formulate the roots of the equation (32) 
as purely imaginary and obtain, by setting the real and the imaginary 
constituent of the equation equal to zero, the conditions 


aQ - —• aQ^ao^ + a^'ca^log id + aQ'^log^co + bQ 1 ^ o>^ - 

^ c o " ip c o '') a;>2 ” c o ,U)2log 00 - V^f + 03 l 0 s 2a3 ) 


+ e^, = 0 


> ( 35 ) 


a^' ^ 0)^ + a^^rtcD-Hog cd - bgu^ - bQ'u^log cd - 


Cq' ^ cd + d 0 ' log co(l - jod) + = 0 


From them cd would have to be eliminated in every special case 
whereby a conditional equation for any of the parameters s, a, r^, 
q2, cdq, Pq or for the flight velocity v results. As an approxi- 
mation, it is sufficient to investigate the equation (33)* Since the 
approximation (3U) shows in the case of instability, a small additional 
damping, the criteria for instability derived from the following equa- 
tions are necessary but not always sufficient; they are therefore valid 
only with the reservation of checking with the approximation (3^) . When 
they are satisfied, the wing may at any rate be regarded as critically 
supported. Before I set down the most general criterion for instability, 

I shall mention one which is very simple but suffices only for cases of 
great instability. This is the condition that the coefficient Cq becomes 
negative 


v2 


c (5 + P 0 V ) 

pn(2s - p 0 ) 


( 36 ) 
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Written as a function of s, the same condition reads 


s 2 o^ 2 - 2s (p Q + an^ 2 ) + p Q 2 + o^ 2 (a 2 + r 2 + q 2 + p Q v) < 0 
Further limits are yielded by the Routh condition 


or 


<*0 

3 0 


< a 0 — + e - ~ 


'°*0 




(&+ P Q v) 


+ P r 


- 2 P 0 S 


r 


< r‘ + 


p 0 € 


M37) 



(*\ 2 - 2p 0 a ) 


11 + p o(f ■ z ) 


X 


In order to satisfy this equation, it will above all be required 
that v be sufficiently large, that is, sufficiently small. Fur- 

thermore, s must be positive; therewith s - a is positive, too, 
because, as mentioned before, a > 0 would for sufficiently large v 
always result in great instability. One can readily understand that 
this must be so. The centroidal axis of the wing as axis of inertia 
generally lags behind the elastic axis as line of application of the 
directional force. If the centroidal axis therefore lies, in flight 
direction, behind the elastic axis (s - a > 0) , the wing has in its 
upward motion, on the average, a positive angle of attack; the opposite 
is true for the downward motion. Thus the air forces always take effect 
in the direction of the motion and amplify the oscillation; whereas in 
the case s - a <0 the opposite, that is, damping occurs. 

It can easily be confirmed that both degrees of freedom must act 
together for achievement of the oscillation. The calculation always 
yields roots with negative real parts if one of the degrees of freedom 
is suppressed (corresponding to = 00 or q2 = 00 ) . This fact is 

confirmed by the failure of tests undertaken formerly in Gbttinger with 
a wing with only one degree of freedom. 
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The power L produced or, respectively, absorbed by the wing is 
according to previous formulas 


L f = A^pTtv^lf 

l f = Z faa + 2b ' z fa0’ + • • • 



( 38 a) 


If I put temporarily co' = 0 + ico, the mean value of the entire 
wing energy for constant amplitude A is 


E = i mv 2 A 2 Q 

4 


Q = cjc^ 2 + cd 2 + 2b ' (ao^ + saf ) + 


1 ' (ao>Q 2 + sod 2 ) 


(b' 2 + b" 2 ) |^, 2 + q. 2 )oD Q 2 + (r 2 + s 2 )cd 2 


A = Age 7 ^, and hence, in a different form 


-L f = — = i mv 2 AAQ = £ mv3A 2 Q 


dt 2 


l r = - rr- Q 


2p, 


0 


(38b) 


If the oscillation calculation has been carried out, the agreement 
of the values If found by different methods offers therefore a control 
for the calc ulat ion. 

7. Comparison with a test result .- For confirmation of the theory, 
tests with a light wooden wing of 60 cm length and 10 cm chord were 
performed in the small wind tunnel of the Gottinger aerodynamic test 
institute. 
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This wing was suspended on an axis between plane side walls by 
means of springs. The springs acted on cross-shaped metal sheets 
which carried moreover sliding weights for variation of the distance 
from the center of gravity and of the moment of inertia. The suspen- 
sion springs or, respectively, their initial compression transferred 
the directional force or, respectively, the directional moment corre- 
sponding to my formulation to the suspension axis. I shall not discuss 
here further details of the apparatus; I refer instead to the drawing 
figure 7. The purpose of the tests was the determination of the incre- 
ments or decrements of the oscillations; for this it was sufficient to 
plot the oscillation of a point, for instance the suspension axis. 

This was done with the aid of a registering drum kindly put at our 
disposal by the physiological Institute of the University Gottingen. 

A few of the diagrams of damped and increasing oscillations thus 
obtained are reproduced in figure 8. From these diagrams the ratios 

jut 

c = e 10 of amplitudes succeeding one another were determined and were 
compared with the values found by calculation with consideration of the 
frictional damping produced by and a or of the translatory or 

rotary oscillations. Within the considerable limits of error which 
could have been cut down only by a major expenditure of time and means, 
the agreement may be called satisfactory. 

The variation of the numbers was in excellent agreement with the 
theory. An increase or damping of the small oscillations resulted 
according to whether the wing was overweight toward the rear or the 
front; the ratio of amplitudes succeeding one another increased when 
the overweight was increased or when the air was blowing more strongly 
against the wing. A precession of the leading edge of the wing was 
clearly evident when the wing was supported so as to be unstable. This 
corresponds in the usual notation to a lag of the rotary oscillation 
as must be the case for power absorption. 

Finally, it is noteworthy that for small air velocities of about 
5 m per second the influence of the viscosity was shown by the fact that 
the lift was not yet fully developed at small angles of attack . For the 
smaller characteristic values of about 500 m/sec. mm a symmetrical wing 
has a small range "dead angle of attack" where the lift remains near 
zero (as shown also by other tests in Gottingen) . The small oscillations 
were therefore still damped when oscillations with a larger initial 
amplitude were already increasing. The limiting amplitude lying between 
these two cases decreased of course more and more with increasing v. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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TABLE I 

NUMERICAL TABLE OF CALCULATED COEFFICIENTS 


0) = 

0.02 

0.04 

0.06 

0.08 

0.10 

0.12 

v = 

0.00262 

0.00773 

0.01351 

0.01892 

0.02317 

0.02572 

*a" = 

0.03852 

0.07387 

0.10605 

O.13532 

0.16182 

0.18601 

y = 

O.O386O 

0.07449 

0.10806 

0.13982 

0.17029 

0.19989 

ky" = 

-0.00184 

-0.00468 

-0.00694 

-O.OO788 

-0.00710 

-0.00452 

= 

-0.0002 

-0.0008 

-0.0018 

-O.OO32 

-0.0050 

-0.0072 

^ It 

0 

0 

0 

0 

0 

0 

niy ' = 

-0.000003 

-0.000024 

-0.000081 

-0.000192 

-0.000375 

-0.000648 

my" = 

0.0004 

0.0016 

O.OO36 

0.0064 

0.0100 

0.0144 

“"^cax = 

-0.9336 

-0.8684 

-0.8073 

-0.7514 

-0.7014 

-0.6576 

’ = 

0.02369 

0.03209 

0.03480 

0.03419 

0.03150 

0.02752 

“ -2w a 7 " = 

-0.4529 

-0.4094 

-O.3707 

-0.3371 

-O.3089 

-0.2862 

OT^Vyy = 

0.02924 

0.05414 

0.07468 

0.09088 

0.10273 

0.11024 

00-2 l cta = 

0.02944 

0.05494 

0.07648 

0.09408 

0.10773 

0.11744 

a) " 2l a7’ = 

O.OOO63 

0.00251 

0.00566 

0.01005 

0.01571 

0.02262 

CD-2l a-y" = 

1 






a/ 

yo. 02944 

0.05494 

0.07648 

0.09408 

0.10773 

0.11744 

ar^Z^ = 

J 






“"^faa = 

0.9631 

0.9233 

0.8837 

0.8455 

0.8091 

0.7750 

co ” 2z fa7' = 

-0.02306 

-0.02958 

-0.02914 

-0.02414 

-0.01579 

-0.00490 

cu ‘ 2l fa7 n = 

0.48235 

0.4643 

0.4471 

0.4312 

0.4166 

0.4036 

O 

03" Ifyy = 

0.0002 

0.0008 

0.0018 

0.0032 

0.0050 

0.0072 
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Figure 2.- Vectorial lift and moment coefficients for equal beat amplitudes. 
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Figure 3.- Vectorial lift and moment coefficients for equal beat velocity. 
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Figure 5.- Degrees of efficiency as a function of co for different amplitude 

ratios c. 
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Figure 6.- Degrees of efficiency as a function of the phase cp. 
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Figure 7.- Test arrangement. 
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Figure 8.- Oscillation diagrams. 
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